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ABSTRACT 

Effect of impact ionization on propagation of a millimeter (mm) 
wave with a Gaussian profile is studied under paraxial ray 
approximation. It is found that early in time the charge density is 
less and hence the defocusing of mm wave is less however later 
in time and space as the carrier density builds up due to impact 
ionization the defocusing of the mm wave is rapid and this 
results in decrease in the intensity of the wave. 
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Interaction of electromagnetic wave 
with semiconductor is an active area of 
research for almost last four decades and 
several nonlinear effects viz. harmonic 
generation, instabilities, wave mixing, etc. 
have been observed 1 ' 2 . When an 
electromagnetic wave passes through a 
semiconductor the valence electron can 
jump to conduction band at the expense of 
the wave energy via one of the following 
processes: impact ionization, tunnel 
ionization or avalanche effect 3 . Some 
notable applications of this transition are in 
switches, millimeter wave generation, UV 
lasers etc. 4 ' 5 . In impact ionization an electron 



or hole can gain energy in the presence of 
electric field and subsequently lose their 
energy by creation of other charge carriers. 
This can lead to avalanche breakdown in 
semiconductors. The electron hole (e-h) 
plasma created via impact ionization is space 
time varying and alters the propagation 
dynamics of the electromagnetic wave. 

In such processes the frequency 
range of incident radiation is important as 
transit time effect come into play. Amongst 
several semiconductor materials, Gallium 
Arsenide (GaAs), a direct band gap 
semiconductor, has some important 
applications such as in microwave frequency 
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integrated circuits, infrared light emitting 
diodes, laser diodes, mobile phones, satellite 
communications, microwave point to point 
links, radar systems and solar cells 6 ' 7 . In the 
microwave frequency range GaAs is one of 
the best candidates because of its higher 
saturated electron velocity, higher electron 
mobility, less noise as compared to Silicon 
(Si) and also can be operated at higher 
powers due to higher breakdown voltages 8 ' 9 . 
In this communication we develop a theory 
to study the effect of impact ionization on 
propagation of a Gaussian millimeter wave 
propagating through GaAs. The wave 
equation is solved by generalizing the 
paraxial ray theory of nonlinear wave 
propagation. 

Consider the propagation of a 
millimeter wave with Gaussian intensity 
profile through a GaAs semiconductor with 
<100> orientation, 

E = E Q (z,r,t)e- i<m - kz) (1) 
Atz = 0 



E Q =A Q e~ r2 ' r \ forO<t<T 



=0 



, otherwise 



(2) 



where x is the pulse duration. 
The wave creates electron - hole (e - h ) 
pairs inside semiconductor via impact 
ionization. The ionization rates for electrons 
and holes can be written as 



dt 
and 



-=a : n 



(3) 



dn 



dt 



i hO ' 



(4) 



respectively. Here n e(h) is electron (hole) 
density, n e(h )o being the equilibrium density, 
and CC i (J3 i ) is electron (hole) ionization rate 
given by 



a^a^e 



E 0n /E 



-E 0p IE 



(5) 
(6) 



a 0 , j3 0 , E 0n and E 0p are constants (c.f. 
Shur 9 pp.188). 

The electron (hole) current density J le(h) is 
governed by 



dJ 



le(h) 



n 



dt 



+ VJ 



'e{h)e 2 E 



le(h) 



(7) 



m 



e(h) 



where v is the electron - hole collision 
frequency and m* e(h) is the electron (hole) 
effective mass. 

The millimeter wave field in the space - 
time evolving electron - hole plasma could 
be written as 

E = Ae*, (8) 

where A is a slowly varying function of z,t 
and <j) is a fast varying function of z, t. 

The wave equation governing the 
propagation of millimeter wave is written 



as 



10, n 



2 - £ L dE _4xdJ 



V l E 



c 2 dt 2 c 2 dt 



(9) 



where e L is lattice permittivity and we have 
neglected the v(v.i?) term for transverse 

waves. For E = yE y form Eq.(8) we have 
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d 2 E 



dz 2 



dz 1 



Ae l * + 



n .d(f>dA 

2i———e 
dz dz 



(10) 



with / as beam width parameter. Using 
Eq.(15) we expand Eqs.(3) and (4) around 
r=0 and obtain 



dm 



-=«o<o 



dt 2 



dt 2 



dj)_ 

K dtj 



Ae* + 



dt dt 
V 2 ± E=V 2 ± Ae», 



(11) 



(12) 



dt 
and 



da? h 



Q e (r = 0) + 



dr 2 



,(16) 



dt 



Q h (r = 0)+ 



dr 2 



,(17) 



respectively. 
Here 

Q e(h) =exp[-(E Qn(p) /E m )(fe r2,2r ° 2f2 )] 



where we have neglected d 2 A/dz 2 terms. 
Defining C0,kas G)=-d0/dz and On integrating Eqn.(16) and (17) we get 

k=d<pldt, with G) 2 = G) 2 p (r = 0) + k 2 c 2 , 

and using Eqs.(7), (11) & (12) in Eq.(9) we 
get 



dA 



+ V 2 L A + 



licoe L dA 



3fc 

~dz~ 



,£ L d(0 



zlv. 



dz ± c 2 dt 

~ ^ l< + < - {col + < L 0 V = 0. d3) 

* respectively 
where (0 2 pe(h) =^Jtn e(h) lm e(h) is the electron Here 

(hole) plasma frequency and m* e , h) is the 

electron (hole) effective mass. 
Multiplying Eq.(13) by A, we obtain 



<= \(a Q -a 0 r 2 /r 2 )dt\ 

zlv g 

and 

t+z/v g 

«& = p Q -0,r 2 lr 2 ]df, 



(18) 



(19) 



■ _ i^ Qn . _ 2 -(E a JE m )f 



dA 2 dA 2 ic 2 



dt 



9z « 26; dt 



a' F 



\a 2 P e + a 2 P h ) 



,2 



r=0 



0.(14) 



Using Eqs.(18) & (19) in Eq.(14) we get 



Later we will recast Eq.(14) in terms of new 
variables z\r,f where t'=t-z/v g ,z'= Z. 

We assume a Gaussian ansatz for the r- 
profile of laser intensity 12 , 

7 2 



dA 



(a p +/3 p ), 



2ik — +V'A+ 

dz' c 2 

Introducing an eikonal 



-A=0 



(20) 



\e \ 2 =— °° * rl: "- 



r 



(15) 



A — Ag(r,z)exp[— S(r,z)] 13 and separating 

the real and imaginary parts of above 
equation we obtain 
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2 — +- 

dz' k 



dr 



(21) 



r 0 k c 



and 

dz' k k dr dr 

respectively. 

Following Eq.(15) we write 
^=^exp(-r 2 /r 0 2 / 2 ), 
and expand 5 as 

S=|y0(zVW), (24) 

in the paraxial ray approximation. Using 
Eqs.(23) &(24) in Eqs.(21) & (22) we get 
1 df 



(23) 



P= 



fdz' 



(25) 



and the equation governing the beam width 
parameter /as 



1 



-+ 



(a p +j3 p )fr 2 



2 2 



(26) 



dz' 2 R 2 f 
where R 2 d =k r 0 2 . 

We introduce dimensionless variables 
g = z'/R d and 7=10 



2 <^ f then 



2c 2 



Eq.(26) could be rewritten as 
df f 



If 

o J 



2 

■^00 MpeO 



drj 



(27) 



We have solved Eq.(27) numerically for 
E 0n /E m =6.27 ,E 0p /E m =4.S2, 



n e ~n h ~ 5X10 6 cm~ 3 ,e L ~ 11, 

X = 1cm, r Q = Imm, 

m* e = 0.063m 0 ,m* h = 0.063m 0 , 

m 0 = 9.1xl0~ 31 %, tf 0 =2.19xl0 6 , 



(22) P 0 = 2.47x10 6 with the boundary 
conditions: f =1 and df I d% = 0at ^=0 for 
all n. We have chosen time step size 
Ar/ = 0.0 land the space step size 
A£ = 0.01. We write Eq.(27) as 



f"=y +p f 

where 
P = 10 2 . 



(28) 



1 



0 J 



Z7 — — 
a 0 ^On E Ci 



'-f + fi 0 E ( 



0 ^0p °>ph0 E K 

"1 ^ 



CO 



peO 



dri 



and the prime represents differentiation with 
respect to ^. We begin by evaluating P at 
^=0 for all values of T|. We solve Eq.(28) 
with the Runge-Kutta method using the 
value of P obtained at ^=0. Using these 

values of /, we evaluate P at ^ ~ ^ at 
different values of r\. This way we advance 
in \. In Fig. 1 we have plotted the beam 
width parameter / and in Fig.2 the axial 

/ = / / f 1 

intensity 0 J as a function of \ for 
different values of n. Early in time, the axial 
intensity decreases due to diffraction 
divergence. However, as the e-h density 
builds up the effect of self-defocusing 
becomes important and intensity falls off 
more rapidly with the distance of 
propagation. It can be seen that the 
millimeter wave can be severely defocused 
due to impact ionization. 
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Fig. 1. Variation of beam width parameter / with \. 
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